AD-A093  982 


DYNAMICS  TECHNOLOGY  INC  TORRANCE  CA  F/6  1/3 

analysis  for  obtaining  high  lift  pressure  distributions  for  tra— ETC<U> 

JUN  00  E  JAMES r  K  KUSUNOSE  N000U-79-C-0*S8 


Dynamics  Technology,  Inc. 

~ DT-7817-6 


ANALYSIS  FOR  OBTAINING  HIGH 
LIFT  PRESSURE  DISTRIBUTIONS 
FOR  TRANSONIC  AIRFOILS 

Final  Report 


June  1980 

Sponsored  by:  Office  of  Naval  Research 

800  North  Quincy  Street 
Arlington,  Virginia  22217 

Under:  Contract  No.  N00Q14-79-C-0458 

By:  E.  James,  Ph.D. 

K.  Kusunose,  Ph.D. 


Dynamics  Technology,  Inc. 

22939  Hawthorne  Blvd,,  Suite  200 
Torrance,  California  90905 
(213)  373-0666 


DTIC 


D 


Approved  for  Public  Release;  Distribution  Unlimited 


This  report  has  undergone  an  extensive  internal 
review  before  publication,  both  for  technical 
and  non-technical  content,  by  the  Program  Manager 
and  an  independent  internal  review  committee. 

Program  Manager: 


Internal  Review: 


A ^  r'  if 

rT:  ’ 

rnr 

Umnn< 
i  Just  11 

.ion  For 

GRA&I  X 

cab  X] 

iiinced  □ 

fty  -  -  -  - 

Distribution/ _ 

Availability  Codes 

Avail  and/or 

DlBt 

Special 

DTlC 

electe 

JAN  2  2  1981 


SECURITY  CLASSIFICATION  OF  This  PAGE  rWhtn  Of  Cntmd) 

I  REPORT  DOCUMENTATION  PAGE 


OPT NUMBCW 

7R17-fi  / 


12.  govt  accession  no.i 


i  PT-7817^6 ,_C-  _ }/?  II- ft, 

4.  TITLE  *n»nr/»>  - 

^Analysis  for  Obtaining  High  Lift  Pressure 
Distributions  for  Transonic  Airfoilsf# 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3.  SccSIcSrs  CATALOG  NUMBER 


5.  tyi »€  or  nt*o*r-±  »cmoo-CQvCRCO 


^  Technical  Report 


•i  / 


17.  author*; 


«.  PCnrORMINO  ORG.  REPORT  NUMBER 

DT-7817-6 

a.  CONTRACTOR  GRANT  NUMflERf *) 


E.yjames,  K.^<usunose 


■79-C-0458 


».  performing  organization  name  ano  aooress 

Dynamics  Technology,  Inc.  j,  M  Y 

22939  Hawthorne  Blvd.,  Suite  200  )  >>  ^ 

Torrance,  California  90505 _ -  -  . _ 

II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research,  Dept,  of  the  Navy  // 
800  North  Quincy  Street  5 

Arlington,  Virginia  22217 

<4.  MONITORING  AGENCY  NAME  A  ADDRESS///  dlllHn I  /ran  CantraJI/n.  OfficaJ 


10.  program  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 


NR  212-263 


12.  REPORT  OATE 

June,  198^ 

IT  NUMflER  OP  PACCfi  - 
91 

15.  SECURITY  CLASS,  fair  rhia  roportf 


UNCLASSIFIED 


1S«.  OeCLASSlPlCATlON/  oowmgaaoing 
SCHEDULE 


[  16.  OlSTPlBUTJOM  STATEMENT  foi  thia  Raport) 


Approved  for  public  release;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (ol  thm  abatract  at uarad  in  Block  30,  It  dlttarant  Irom  Report) 


\  15.  SUPPLEMENT A*Y  NOTES 


19.  KEY  WORDS  fContinua  on  ravaraa  at  da  it  nacaatary  and  Idanlltv  by  block  numbar) 

transonic,  airfoil,  design,  maximum  lift,  pressure  distribution, 
compressible  separation  criterion,  shocks,  variational  problem, 
turbulent  boundary  layer 

20  UBSTRACT  [Continue  nn  r*w»  P4#  If  n.c.ny  IRnflfv  6v  block  number) 

vVhis  report  describes  the  development  of  a  method  for  obtaining  an  optimum 
pressure  distribution  over  the  suction  side  of  an  airfoil  operating  in  the 
transonic  speed  regime.  The  pressure  distribution  is  optimum  in  the  sense 
that  maximum  lift  on  the  airfoil  section  is  desired  while  the  flow  over  the 
airfoil  remains  fully  attached. 

The  development  allows  for  shocked  and  shock  free  flows  over  the  airfoil  with 
a  compressible  and  mixed  laminar-turbulent  boundary  layer.  The  boundary  layer 


security  Classification  OF  this  page  .■**•«  £>•/•  £ni» rad. 


DO  ,  1473  EOlTlON  OF  I  NOV  SS  IS  OBSOLETE 


umrv  classification  of  this  FAoemM*  Dim  bim« 


is  typically  maintained  on  the  verge  of  incipient  separation  over  the  rear¬ 
ward  stretch  of  the  airfoil.  The  pressure  recovery  distribution  is  obtained 
from  a  compressible  extension  of  Stratford's  incompressible  turbulent  boundary 
layer  separation  criterion  which  is  also  given  in  the  report. 

A  variational  problem  for  maximum  lift  is  formulated  and  solved  to  yield  the 
location  of  the  onset  of  the  pressure  recovery  and  the  pressure  distribution 
forward  of  this  point.  _ 


HCUWTY  CLASSIFICATION  OF  This  FAaer»Ti»n  0««  Cm.no 


PREFACE 


This  work  was  supported  by  the  U.  S.  Office  of  Naval  Research  under 
Contract  No.N00014-79-C-0458.  The  authors  would  like  to  express  their 
appreciation  to  Or.  A.  M.  0.  Smith  (Consultant)  and  Dr.  R.  Whitehead  (of 
ONR)  who  helped  shape  our  thoughts  on  the  subject. 


ABSTRACT 


This  report  describes  the  development  of  a  method  for  obtaining  an  opti¬ 
mum  pressure  distribution  over  the  suction  side  of  an  airfoil  operating 
in  the  transonic  speed  regime.  The  pressure  distribution  is  optimum  in 
the  sense  that  maximum  lift  on  the  airfoil  section  Is  desired  while  the 
flow  over  the  airfoil  remains  fully  attached. 

The  development  allows  for  shocked  and  shock  free  flows  over  the  airfoil 
with  a  compressible  and  mixed  laminar-turbulent  boundary  layer.  The 
boundary  layer  is  typically  maintained  on  the  verge  of  incipient  separa¬ 
tion  over  the  rearward  stretch  of  the  airfoil.  The  pressure  recovery 
distribution  is  obtained  from  a  compressible  extension  of  Stratford's 
Incompressible  turbulent  boundary  layer  separation  criterion  which  is 
also  given  in  the  report. 

A  variational  problem  for  maximum  lift  Is  formulated  and  solved  to  yield 
the  location  of  the  onset  of  the  pressure  recovery  and  the  pressure  dis¬ 
tribution  forward  of  this  point. 
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comparison  flow  (usually  taken  as  6  or  7  for  shock 
free  flows) 

local  pressure 

pressure  in  the  recovery  region  along  the  edge  of  the 
boundary  layer  downstream  of  the  reference  point 

reference  pressure 

nondlmenslonal  pertubation  pressure 

Reynolds  number  based  on  u0,  vQ  and  s;  uQs/vo 

Reynolds  number  based  on  uj_,  n^  and  s;  u^s/v^ 

Reynolds  number  based  on  u0,  and  s0;  U0SQ/V0 

Reynolds  number  based  on  uj.nj  and  ss0;  u^sSQ/v1 

Reynolds  number  based  on  ;u0,  vQ,  and  x0;  uQx0/v0 

free  stream  Reynolds  number  for  unit  chord  length  s  U 

recovery  factor  (usually  =*0.9) 
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*  n-2 
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the  shock  case  p  1 
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equivalent  shock  free  flat  plate  length  for  the  shock  case 
equivalent  flat  plate  trailing  edge  length 
absolute  temperature 

velocity  component  In  free  stream  direction 

comparison  flat  plate  velocity  profile,  uc( s , 
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u  at  the  trailing  edge,  ue(l) 

defined  by  (21) 

transverse  velocity  component 

arc-length  measured  along  the  top  airfoil  surface 
starting  at  the  nose 

x  at  the  maximum  velocity  point  on  the  airfoil 

local  measure  of  transverse  distance  along  the  surface 

y  at  the  join  of  inner  field  and  far  field 

defined  by  (A-3) 
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boundary  layer  thickness 
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constant  defined  by  (A-10) 

speed  ratio  scaled  with  the  reference  speed  u0,  u/u0 
defined  by  (A-22),  uc(s,1>)/u0 
constant  defined  by  (A-10) 
momentum  thickness 

shock  Induced  momentum  thickness  defined  by 


von  Karman  constant,  0.41 

pressure  gradient  coefficient  |p^k:~2u~2 

coefficient  defined  by  the  power  law  expression  for  the 
comparison  flow  (A-30) 

defined  by  (B-3) 

kinematic  viscosity,  u/p 

constant  defined  by  (A-10) 

fluid  density 

fluid  density  along  the  edge  of  the  boundary  layer 
fluid  density  along  the  wall 
defined  by  s/s0 

nondlmensional  perturbatin  density 
shear  stress 
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total  velocity  potential 
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NOMENCLATURE  ( Conti nued) 


$  nondimensional  perturbation  velocity  potential 

stream  function  for  compressible  flow  defined  by  { A- 1 3 ) 

^  limiting  streamline  on  the  Inner  periphery  of  the  outer 

regi  on 

(  )a  refers  to  adiabatic  conditions 

(  )c  refers  to  conditions  of  comparison  flat  plate  flow 

(  )e  refers  to  conditions  along  the  edge  of  the  boundary  layer 

(  )tr  refers  to  conditions  at  the  transition  point 

(  )w  refers  to  conditions  along  the  wall 

(  )0  refers  to  reference  conditions;  also  refers  to  conditions 

ahead  of  the  shock  for  the  shock  cases 

(  refers  to  conditions  aft  of  the  shock 

(  )„  refers  to  freestream  conditions 


1.  INTRODUCTION 


For  an  aircraft  to  maneuver  effectively  in  the  transonic  speed  regime, 
its  wing  must  develop  high  lift  without  incurring  excessive  drag.  This 
report  presents  a  methodology  for  developing  high  lift  on  airfoil  sec¬ 
tions  by  seeking  a  pressure  distribution  that  achieves  high  lift  while 
maintaining  fully  attached  flow  along  the  airfoil  surface.  Conse¬ 
quently,  the  importance  of  reliably  predicting  incipient  separation  over 
an  airfoil  section  operating  In  the  transonic  speed  regime  and  In  iden¬ 
tifying  and  understanding  how  the  principal  controlling  parameters  of 
foil  shape  and  flow  characteristics  affect  separation  have  been  identi¬ 
fied  as  a  principal  area  of  this  investigation.  This  Issue  is  particu¬ 
larly  relevant  to  the  primary  objective  of  developing  a  rational  method¬ 
ology  for  designing  a  high  lift  wing  section  for  transonic  maneuver¬ 
ability.  Activities  have  therefore  concentrated  upon  the  analysis  of  a 
turbulent  boundary  layer  and  its  separation  characteristics  over  the  low 
pressure  side  of  a  transonic  wing  section.  In  this  endeavor,  a  compres¬ 
sible  separation  criterion  has  been  developed  In  which  heat  transfer 
effects  have  also  been  Included.  The  criterion  Is  similar  to  that  of 
Stratford's  Incompressible  turbulent  separation  criterion  [I], 

The  procedure  for  determining  the  high  lift  pressure  distribution  fol¬ 
lows  closely  that  of  Liebeck  and  Smith  [2]  which  has  led  to  successful 
high  lift  incompressible  airfoil  shapes.  In  fact,  the  variational 
problem  for  maximum  lift,  results  in  all  cases  (incompressible,  compres¬ 
sible,  with  and  without  shocks)  to  the  roof-top  pressure  distribution 
followed  by  a,  so  called,  Stratford  type  pressure  recovery.  The  loca¬ 
tion  of  the  on-set  of  the  pressure  recovery  and  the  shape  of  the 
recovery  pressure  distribution  depend  upon  the  free  stream  conditions. 

Maximum  lift  is  achieved  by  the  shock  free  case  for  specified  free 
stream  conditions.  For  the  same  free  stream  conditions,  a  shock  may 
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occur  which  will  result  In  a  reduction  of  lift  (compared  with  the  shock 
free  optimum  pressure  distribution)  and  an  increase  in  drag. 

An  effort  Is  anticipated  In  the  near  future  to  Incorporate  the  pressure 
distributions  obtained  herein  Into  a  transonic  airfoil  design  code  to 
generate  high  lift  airfoil  shapes.  That  effort  will  concentrate  on 
obtaining  practical  shapes  from  which  a  high  lift  maneuvering  wing  can 
be  designed. 

Inasmuch  as  maximum  lift  also  depends  on  the  specification  of  a  realis¬ 
tic  pressure  distribution  on  the  lower  surface  of  the  airfoil,  the 
follow-on  effort  to  this  work  will  Initially  concentrate  on  the  develop¬ 
ment  of  such.  Determination  of  an  appropriate  lower  surface  pressure 
distribution  when  the  upper  surface  distribution  is  specified  brings  up 
the  matter  of  finding  a  compatible  pressure  distribution  for  transonic 
airfoil  design  codes.  That  Is,  In  determining  one  that,  with  given  free 
stream  conditions,  will  satisfy  the  closure  condition  at  the  airfoil's 
trailing  edge. 


2.  COMPRESSIBLE  TURBULENT  BOUNDARY  LAYER  SEPARATION  CRITERION 


A  compressible  turbulent  boundary  layer  separation  criterion  is  sought 
using  the  line  of  reasoning  of  Stratford  in  his  original  formulation  [1] 
of  an  incompressible  turbulent  boundary  layer  separation  criterion.  His 
criterion  is  sufficiently  accurate  and  Is  convenient  for  practical 
applications.  These  traits  motivated  the  present  attempt  for  an  exten¬ 
sion  to  account  for  compressibility. 

The  extension  appears  possible  giving  rise  to  a  criterion  having  Mach 

number  dependence  on  the  right  hand  side  of  the  external  pressure  and 

longitudinal  pressure  gradient  expression.  The  form  is  similar  to 

Stratford's  but  with  slight  changes  in  certain  empirical  constants,  and 

in  the  definition  of  an  appropriate  pressure  coefficient  for  a  compres- 
★ 

sible  fluid,  Cp.  This  latter  quantity  is  defined  as 


where  the  subscript  'o'  denotes  a  reference  condition.  For  an  airfoil 
surface,  the  reference  condition  corresponds  to  the  point  on  the  surface 
where  the  Inviscld  solution  achieves  minimum  pressure  (maximum  stream- 
wise  flow  speed). 

y  is  the  ratio  of  specific  heats  of  the  fluid. 

p0  is  the  reference  pressure. 

u0  Is  the  streamwise  flow  speed  at  the  reference  point. 

pe  denotes  the  pressure  In  the  recovery  region  along  the  edge  of 
the  boundary  layer  downstream  of  the  reference  point.  This 
pressure  is  Impressed  on  the  boundary  layer  and,  therefore,  is 
assumed  not  to  vary  with  transverse  position  within  the 
boundary  layer. 

ue  Is  the  streamwise  flow  speed  along  the  edge  of  the  boundary 
layer. 
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The  turbulent  compressible  Stratford  type  separation  criterion  was  found 


to  be 


/.  < 


(•*)  ->ft)  ( 


R  .10 
s 


where 


g(n)  =  <  k1n/2b(n) 


b(n)  s  21.32 


n-2 

Tt  3  \1/4 


n2(n+l) ■ 


s  Is  arc-length  measured  along  the  airfoil  surface  based  upon  an 
equivalent  flat  plate  length.  To  obtain  the  origin  of  s  the 
following  procedure  is  used: 

(I)  Locate  at  the  edge  of  the  boundary  layer  the  maximum 
flow  speed  uQ. 

(II)  Determine  the  boundary  layer  momentum  thickness  at  that 
point. 

(ill)  Consider  flow  past  a  semi -infinite  flat  plate  at  zero 
incidence  with  free  stream  conditions  corresponding  to 
Pq ,  u0,  Pq .  Determine  the  length  of  flat  plate  Sg 
necessary  to  duplicate  the  momentim  thickness  achieved 
by  the  airfoil  where  u0  occurs. 

(iv)  The  origin  of  s  is  taken  to  be  the  equivalent  flat 
plate  distance  s0  upstream  from  the  point  of  maximum 
flow  speed  on  the  airfoil. 

n  is  the  inverse  of  the  exponent  of  the  power  law  assumed  for  the 
comparison  flow  (usually  taken  as  6  or  7  for  shock  free  flows). 

<  is  the  von  Kinnan  constant  =  0.41. 

k1  is  an  empirical  constant  obtained  from  turbulent  boundary  layer 
flat  plate  experiments  =  0.684. 

pw  is  the  fluid  density  along  the  wall. 

Pe  is  the  fluid  density  along  the  edge  of  the  boundary  layer. 

Rs  is  the  Reynolds  number,  based  on  u0,  s  and  (the  kinematic 
viscosity  of  the  fluid  at  the  reference  point). 


The  density  ratio  (assuming  the  wall  Is  adiabatic)  IpJpA  introduces  a 

e  a 

Mach  number  dependence  Into  the  separation  criterion  since 

3 1  +  r  Me  (3) 

where 

r  is  the  recovery  factor  (usually  •  0.9),  and  for  turbulent 
boundary  layers  It  Is  the  cube  root  of  the  Prandtl  number. 

Me  Is  the  Mach  number  based  on  ue  (Mach  number  based  on  the  Invls- 
cld  solution,  and  therefore  dependent  on  location  s). 

Appendix  A  of  the  report  presents  a  detailed  analytical  development 
leading  to  the  separation  criterion  just  given.  Use  Is  made  there  of  an 
outer  boundary  layer  flow  corresponding  to  a  region  near  the  edge  of  the 
boundary  layer  and  an  inner  flow  region  very  near  the  surface  of  the 
airfoil.  Representations  of  the  flow  In  these  two  limiting  boundary 
layer  regions  are  then  patched  together  In  an  Intermediate  region  to 
arrive  at  the  separation  criterion. 

2.1  Procedural  Stannary  for  Obtaining  the  Separation  Criterion 

In  the  derivation  of  the  separation  criterion,  the  Prandtl  mixing  length 
hypothesis  Is  used  to  represent  the  shear  stress  in  the  inner  region 
close  to  the  wing  surface.  An  expression  for  the  shear  stress  Is  also 
provided  by  Integrating  the  boundary  layer  equations  outward  from  the 
surface  to  the  point  y  of  the  Inner  region  where  the  stress  Is  to  be 
evaluated.  Equating  the  two  expressions  gives  a  nonlinear  first  order 
ordinary  differential  equation  for  the  velocity  profile  In  the  Inner 
region.  The  solution  (which  can  be  obtained  in  closed  form)  has  the 
following  asymptotic  form  valid  close  to  the  surface  when  the  wall  shear 
stress  vanishes: 
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S-.  (l  *  $  *  M)  ; 

where  .  ,, 

A  -  £  £  )1/2  (A. 

81  *  (l  *  ^  ^)fe)-  1 

Me  Is  the  Mach  number  at  the  edge  of  the  boundary  layer. 

Te  Is  the  temperature  of  the  fluid  at  the  edge  of  the  boundary 

layer. 

Tw  Is  the  temperature  of  the  airfoil  surface. 

Thus,  when  the  turbulent  boundary  layer  Is  on  the  verge  of  separating, 
its  velocity  profile  tends  to  zero  as  the  square  root  of  the  distance 
from  the  surface.  Equation  (4)  characterizes  the  turbulent  boundary 
layer  in  the  inner  region  when  the  flow  is  about  to  separate. 

In  the  outer  region  of  the  turbulent  boundary  layer,  it  is  assumed  that 
the  total  pressure  loss  along  a  streamline  is  independent  of  the  pres¬ 
sure  rise,  a  result  supported  by  experimental  results  (see  Grabowski , 
@t  al.  [3]).  As  a  consequence,  the  pressure  loss  along  a  stream  line  in 
the  actual  case  is  very  nearly  the  same  as  on  a  corresponding  streamline 
in  the  turbulent  flat  plate  case  where  the  pressure  is  constant.  The 
characterization  of  the  velocity  profile  in  the  outer  region  is  given  by 
the  Bernoulli  equation  applied  along  streamlines  where  the  pressure  is 
assumed  related  to  the  density  isentropically. 


To  evaluate  the  constant  of  the  Bernoulli  equation,  a  comparison  flat 
plate  velocity  profile  uc  Is  assumed  to  exist  (similar  to  the  actual 
profile)  and  such  that,  at  the  location  along  the  airfoil  surface  cor¬ 
responding  to  the  beginning  of  the  pressure  recovery  region,  p  3  p0  and 
u  *  uc.  At  points  downstream  where  p  >_  p0,  a  power  law  expression  Is 
assumed  for  uc.  This  completes  the  description  of  the  outer  region  of 
the  turbulent  compressible  boundary  layer  on  the  verge  of  separation. 

The  Inner  and  outer  representations  of  the  velocity  profile  join  at  some 

3C  3 

Intermediate  location  yj  determined  by  calculating  the  quantity  'K-gpO 
In  each  of  the  regions  and  equating.  Here,  is  the  streamf unction  and 
C=u/u0  Is  the  speed  ratio.  The  relationship  between  u  and  uc  at  the 

join  is  determined  by  a  similar  procedure  from  the  expression 
.  The  separation  criterion  (2)  then  follows  from  the 
Bernoulli  equation. 
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3.  COMPRESSIBLE  STRATFORD  FLOWS 


In  the  previous  section,  a  compressible  separation  criterion  was  given 

for  turbulent  boundary  layers  (equation  (2)).  From  equation  (2)  It  is 

possible  to  derive  theoretical  boundary-layer  flows  downstream  of  the 

peak  velocity  point  where  the  flow  is  on  the  brink  of  separation  or  free 

of  surface  shearing  stress.  Such  flows  are  obtainable,  initially,  by 

★ 

integrating  the  separation  criterion,  from  s  *  s0  where  Cp*  0  to  s  *  sc 
where  Cp3  (n-2)/(n+l).  At  s  *  sc  the  inner  region  of  the  boundary 
layer  reaches  all  the  way  to  the  outer  edge  of  the  boundary  layer.  For 
s  >_  sc,  a  momentun  integral  form  is  used  to  obtain  the  Stratford  flow 
which  continues  the  one  obtained  from  the  separation  criterion. 


Appendix  B  presents  a  detailed  derivation  of  the  Stratford  flows  and  the 
momentum  integral  form  used  to  continue  the  flows  beyond  s  *  sc. 


Treating  equation  (2)  as  an  ordinary  differential  equation  and  special¬ 
izing  to  the  case  when  n  *  6  (i.e.,  the  velocity  profile  follows  a  l/6th 
power  law),  we  obtain 


*  r  /  s  \1/6  i1/3 

cp(s)-[b(|-)  -ij  ;  s<sc 

Cp(s)  -  1  -  (v  +  a4r1/2  ;  S  >  sc 


I 

i 


(5) 


where 


B  =»  1.695  (l  +  r  1/3 (r$  •  10'6) 


.1/18 


R  3  u  s  /v 
s  0  0  0 

*  ■  [>  *  w3]6 


-  1/6  r  1/6  -I 

‘•412W  M  -1. 


■2/3 
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a4  *  5.444  - 


Mq  Is  the  peak  Mach  number  which  occurs  at  s  *  s0. 


■jf  j 

The  quantity  Cp*  l-(ue/uQ)  can  be  related  to  Stratford's  canonical 
pressure  coefficient  (the  pressure  difference  with  respect  to  the  mlnl- 
mun  pressure  p0  normalized  with  respect  to  the  dynamic  head  at  the  mini¬ 
mum  pressure  point,  s  a  sQ).  That  Is, 

rp  -  (P-P0)  /  7  pouo 


The  relationship  (derived  from  equation  (A-48))  is  given  by 


★ 

When  MQ  *  0,  It  can  be  shown  that  Cp  reduces  to  TTp.  That  is, 

C*  ■  IT  ■  1  -  \yZJul  ;  M  =  0 
p  p  e  o  o 


The  compressible  Stratford  flows  (velocity  profiles  on  the  verge  of 
separation)  are  obtained  from  equation  (5)  since 


ue(s) 


(7) 


Figure  1  displays  lTp  versus  C*  (as  given  in  equation  (6))  for  a  range  of 

peak  Mach  numbers  NL.  It  can  be  seen  that  as  the  Mach  number  increases 
V  * 

for  a  given  speed  ratio  (given  Cp)  then  so  also  does  17p.  Further¬ 
more,  l?p  can  exceed  unity  in  contrast  to  the  incompressible  case. 


Figure  2  shows  a  range  of  compressible  Stratford  flows.  As  the  peak 
Mach  number  increases,  compressibility  permits  the  attainment  of  larger 


speed  ratios  (relative  to  the  Incompressible  Stratford  flow)  at  any  par¬ 
ticular  point  s/s0.  Consequently,  for  incipient  separation,  compres¬ 
sibility  allows  the  flow  to  decelerate  at  a  lesser  rate  than  its  incom¬ 
pressible  analogue. 

★ 

Figures  3  and  4  present  the  compressible  Cp  and  adiabatic  density 
ratio  (p/pJ  variations  corresponding  to  the  Stratford  flows  given  in 

a  * 

Figure  2.  It  can  be  seen  that  compressibility  provides  Cp  a  slower 

recovery  than  the  incompressible  case.  Regarding  the  adiabatic  density 

ratio,  pw  decreases  as  the  edge  Mach  number  increases  for  a  specified  p6 

due  to  viscous  heating. 


-15- 


t 


4.  COMPARISON  OF  COMPRESSIBLE  SEPARATION  CRITERIA 

In  addition  to  Stratford' s  original  turbulent  Incompressible  separation 
criterion,  three  formulae  exist  for  the  prediction  of  separation  which 
purport  to  include  the  effect  of  compressibility.  These  formulae  were 
derived  or  suggested  by  Gadd  [4],  Smith  [5]  and  the  criterion  developed 
here  (equation  (2)). 

Assessing  a  number  of  incompressible  separated  flows  showed  that  a 
better  constant  on  the  right  hand  side  of  Stratford's  Incompressible 
criterion  is  0.50  Instead  of  0.39.  At  the  time  Stratford  developed  his 
criterion  there  were  limited  experimental  data  to  aid  in  establishing 
the  constant.  For  purposes  of  comparing  the  different  criteria,  the 
0.50  constant  will  be  used. 


The  relevant  formulae  to  be  compared  are: 
dT_  1/2 


?p  (*  £)  uo-6-«srl/1°  -  o. 

1/2 


50 


Stratford  Incompress. 


Mo  L 


(‘-y-u -a 


*•> 


(l0-6-Rsr1/10=  0.50 


Gadd 


R$r1/10  =  0.50 

„  1/2 

Rsr1/10=0.52  ly-\  (10"6*Rs)'1/6° 


Smith 


(Present) 


(8) 

(9) 

(10) 

(11) 


» 


l 
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These  formulae  have  been  specialized  to  a  particular  exponent  In  a  power 
law-type,  boundary  layer  profile  (n= 6).  In  addition,  Stratford's 
results  have  been  adjusted  to  agree  with  test  data  as  mentioned  previ¬ 
ously.  All  formulae  have  basically  the  same  structure  but  equation  (11) 
contains  a  density  ratio  and  a  slight  Reynolds  number  dependence  on  the 
rlgh  hand  side. 

Since  Stratford's  result  is  strictly  incompressible  and  since  compres¬ 
sibility  affects  Cp  values  (see  Figure  4),  Gadd,  In  view  of  equation 
(9),  proposed  that  compressibility  be  accounted  for  by  the  Mach  number 
ratio  dependence  of  equation  (9).  In  equation  (9),  Me  is  the  edge  Mach 
number  and  M0  is  the  peak  Mach  number. 

Based  on  studies  in  an  AGARD  paper  [8],  Smith  suggested  that  velocities 
be  used  to  eliminate  the  effect  of  compressibility  on  the  pressure  coef¬ 
ficient  and  equation  (10)  Is  the  result.  The  last  equation  (11)  Is  the 
present  version  which  parallels  Stratford's  original  analysis  while 
including  compressibility. 

The  process  for  predicting  separation  by  equations  (8),  (9)  and  (10)  is 
to  observe  the  growth  of  the  left  hand  sides  -  when  the  value  reaches 
0.50,  then  separation  is  assumed  to  occur.  Equation  (11)  is  slightly 
more  complicated  since  its  right  hand  side  is  not  constant.  However, 
the  process  for  predicting  separation  is  essentially  the  same.  A  way  of 
comparing  the  formulae  is  to  consider  various  decelerating  flows  and 
determine  how  the  left  hand  sides  behave.  We  begin  by  converting  equa¬ 
tions  (8)  through  (11)  to  expressions  involving  only  Cp.  For 
Stratford's  formula  we  define 

Cn5  "  c„(M  *° j . 

p  p  p  0 


Using  the  energy  equation,  the  following  relation  can  be  derived  for  use 
in  Gadd1 s  formula: 


1 


From  this  equation,  one  obtains 


s 


d 

Hs 


t1  ^  jr  "o) 

(“T*S  c*p)2 


To  simplify  matters  slightly,  we  introduce  a  quantity  a  *  s/sQ  and  re- 

*  o 

gard  Rs  and  Cp  as  functions  of  o.  Suppose  further  that  R$  *  10  so  that 


u 

_  0  0 
V  — • 


a  x  106 


Consequently,  equations  (8)  through  (11)  can  be  written 


2A 


2A 


2A 


2A 


o 


1 


0 


2 


Stratford  Incompress.  (12) 


Gadd  (13) 


Smith  (14) 


(Present)  (15) 
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where 

Aosl 


Y-l 

~r 


,Y-1 

~T 


1/60  . 
a2  2  -  hr  i1 


+  r 


Y-l 

T 


u 


l1  +JTlMo)(1  +  r  “T‘Mo)”1  "o 


(16) 


To  compare  the  various  separation  criteria,  we  assume  a  linear  velocity 
profile  in  the  form 


—  *  (1  +  6)  -  So 


8  >  0 


The  parameter  8  is  selected  for  this  investigation  to  be  1/3  so  that 
when  a  ■  4,  the  speed  ratio  vanishes.  We  define  a  quantity  G  as 

ar*  x  r  2-,  u  1/2 


G  5  2C_  o 


In  Figure  5,  curves  of  G  and  1/Aj  (1*0, 1,2)  versus  a  are  presented  for  a 

.1 

value  of  Mach  number  M0*1.4  and  8*1/3.  G  intersects  the  curves 
where  separation  is  predicted. 


Gadd's  method  predicts  separation  early  (relative  to  the  other  two 
methods)  for  the  3,M0  constants  selected.  The  Smith  and  present  methods 
predict  separation  very  near  the  same  location.  Consequently,  the  lat¬ 
ter  analysis  tends  to  validate  (and  slightly  refine)  Smith's  original 
contention  for  extending  Stratford' s  criterion  to  the  incompressible 
regime.  On  the  basis  of  the  analysis  developed  herein,  the  present 


James/Smi th  (E«|.  15)  for  a  Linear  Speed  Ratio  of  Slope  -  1/3 
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formula  predicts  separation  to  occur  as  much  as  ten  percent  earlier  than 
the  Smith  formula  when  the  flow  is  such  as  to  Induce  early  separation 
(l.e.,  (s/sQ)  near  unity).  When  s/sQ  >  3,  the  trend  reverses  and 
Smith's  formula  will  predict  a  slightly  earlier  separation  for  those 
flows  such  as  to  induce  separation  for  s/s0  >  3.  The  same  basic  trends 
prevail  with  respect  to  the  Smith  and  Gadd  formulae.  However,  the  per¬ 
cent  variation  is  seen  to  be  considerably  greater. 
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5.  ESTIMATION  OF  LAMINAR  RUN  FOR  COMPRESSIBLE  BOUNDARY  LAYERS 


Lift  generation  on  an  airfoil  section  is  known  to  be  enhanced  by  extend¬ 
ing  the  region  of  laminar  flow  on  the  foil.  The  type  of  flow  we  are 
striving  to  achieve  Is  one  that  rapidly  accelerates  about  the  nose  of 
the  foil  from  the  stagnation  point  x=0  to  some  point  x  on  the  top  sur¬ 
face  (here  x  is  arc-length  measured  along  the  top  airfoil  surface  start¬ 
ing  at  the  front  stagnation  point).  Thereafter  it  only  gradually  accel¬ 
erates  to  its  peak  velocity  which  we  assume  Is  reached  at  x=xQ.  (For 
purposes  of  this  analysis,  It  Is  assumed  that  the  boundary  layer  flow 
remains  laminar  up  to  the  point  xtr  where,  thereafter,  it  spontaneously 
transitions  to  a  turbulent  boundary  layer.)  In  the  region  of  rapid 
acceleration  near  the  stagnation  point  (zone  of  favorable  pressure  gra¬ 
dient)  we  make  the  assumption  that  the  boundary  layer  grows  incompres- 
sibly.  The  momentum  thickness  in  the  laminar  region  can  then  be  calcu¬ 
lated  by  the  method  of  Walz  [6],  which  follows 


,  0.470v 

92 - £ 

ue(x) 


^  /ujC)'5 


,ue(x). 


<15  ; 


x  ixtr 


(17) 


For  the  remainder  of  the  flow  up  to  the  peak  velocity  we  assume 
Ug(x)*0  and  ue*uQ.  This  is  a  flat  plate  assumption  -  but  in  this 
region  we  want  to  account  for  compressibility.  Following  Gruschwitz  [7] 
(see  also  Reference  [8],  Schlichtlng,  pages  341-344),  we  have  for  the 

A 

shape  factor  K*0  corresponding  to  Ug30  that 
■  0.2349 

This  integrates  to 


where  x*  Is  a  constant  of  integration. 


We  observe  that  when  ug(x)  =  u0,  equation  (17)  reduces  essentially  to 
the  previous  equation  when  x*  =  0.  Consequently,  this  result  suggests 
that  we  can  use  equation  (17)  to  a  good  approximation  for  the  growth  of 
a  compressible  boundary  layer  subject  to  a  pressure  gradient  provided 
the  compressible  ue(x)  is  used.  This  same  result  has  also  been  obtained 
from  a  more  rigorous  mathematical  treatment  but  the  analysis  is  lengthy 
and  is  not  provided  here. 

To  determine  the  transition  point  xtr  a  simple  transition  criterion 
based  on  momentun  thickness  Reynolds  number  Re  is  used.  Reference  [9], 
(page  332),  suggests  that  R0  3  400  suitably  determines  transition  for 
compressible  flows  with  pressure  gradient.  When  this  criterion  is  used 
with  equation  (17),  the  laminar  run  can  be  determined. 


5.1  Equivalent  Flat  Plate  Length 

An  equivalent  flat  plate  length  s0  must  be  obtained  for  use  in  the  sepa¬ 
ration  criterion  previously  given.  The  procedure  for  determining  an 
equivalent  flat  plate  length  Is  to  find  the  point  on  a  flat  plate  where 
its  momentun  thickness  Is  equal  to  the  momentun  thickness  at  the  peak 
velocity  point  on  the  airfoil. 

In  this  section,  we  derive  an  expression  for  the  turbulent  momentum 
thickness  which  develops  as  the  boundary  layer  flow  sweeps  o.*er  a  speci¬ 
fied  smooth  surface  from  a  stagnation  point  and  past  the  transition 
point  where  the  flow  is  assumed  to  become  turbulent  Instantaneously.  In 
this  endeavor  we  assume  the  compressible  turbulent  momentum  thickness  8 
can  be  approximated  by  its  incompressible  analogue: 
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where,  for  accelerating  streams,  (see  Reference  [8],  p.  633-634)  n*4, 
d*3.94  and  C  *  0.01475.  In  this  semi -empirical  result,  the  turbulent 
boundary  layer  begins  at  the  point  xtr.  To  obtain  the  constant  C^,  we 
use  the  laminar  momentum  thickness  0tr  expression  derived  by  Vlalz  (i.e., 
equation  (17)).  Equating  the  previous  equation  to  equation  (17)  with  x 
=*  xtr  and  9  *  9tr  gives 


1 


Therefore,  the  momentum  thickness  at  the  maximum  velocity  point  xQ  on 
the  airfoil  is  expressed  as 


The  momentum  thickness  0Q  for  a  turbulent  transonic  boundary  layer  flow 
past  a  flat  plate  has  been  determined  experimentally  [11]  to  be 


9  *  0.022 


(O' 


1/6. 


(19) 


To  determine  the  equivalent  flat  plate  length  s0  of  the  airfoil  (to  the 
point  of  peak  surface  speed  u0) ,  equations  (18)  and  (19)  are  equated  to 
give 
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For  every  point  x  (such  at  0  <_  x  <_  xQ)  there  is  an  equivalent  flat  plate 
distance  s  (such  that  0  s  sQ)  where  the  correspondence  is  estab¬ 
lished  or  requiring  the  momentum  thickness  on  the  airfoil  at  x  be  equal 
to  the  momentum  thickness  on  the  flat  plate  at  s.  In  general,  s  is 
related  to  x  in  a  nonlinear  fashion  in  the  range  0  x  <_  xQ.  However, 
in  the  range  x0^_xj^l,  the  thickness  growth  on  the  airfoil  is 

to  be  constrained  to  be  the  same  as  obtained  from  a  compressible 
Stratford  type  flow  (which  invokes  an  equivalent  flat  plate  analogy  to 
the  airfoil).  Therefore,  in  that  range  s  and  x  are  related  linearly. 
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6.  VARIATIONAL  PROBLEM  FOR  MAXIMA  LIFT 

In  this  section,  the  variational  problem  for  obtaining  the  maximum  lift 
on  a  transonic  airfoil  section  is  formulated  and  solved.  The  theoreti¬ 
cal  pressure  distribution  has  been  obtained  which  provides  the  maximum 
lift  on  the  airfoil  assuming  the  flow  is  fully  attached  over  the  air¬ 
foil.  The  optimum  pressure  distribution  obtained  under  the  constraint 
of  nonseparation  is  comprised  on  the  upper  surface  of  an  initial  flat 
roof-top  pressure  profile  which  prevails  from  the  leading  edge  to  a 

point  xQ  where  the  boundary  layer  flow  is  thereafter  on  the  brink  of 
separation.  Aft  of  this  point,  the  pressure  recovery  is  taken  to  be  a 
compressible  turbulent  Stratford-type  profile.  Such  a  profile  just 
avoids  separation  (within  a  certain  margin  of  safety  to  allow  ;he  flow 
to  remain  attached  under  nominal  off-design  operating  conditions).  On 
the  lower  surface  the  maximum  lift  requirement  imposes  a  stagnation 
(zero  flow)  surface.  Consequently,  the  optimization  problem  has  led  to 
flow  requirements  on  the  upper  and  lower  airfoil  surfaces  that  are 

impossible  to  meet  in  practice.  The  requirements  are  therefore  to  be 
viewed  as  guide  lines  for  achieving  a  high  lift  airfoil  section.  That 
is,  certain  practical  modifications  to  the  theoretically  obtained  pres¬ 
sure  distribution  are  necessary  in  any  realistic  design,  to  achieve 

physically  useful  airfoil  shapes.  These  modifications  are  indicated 

later  in  this  report. 

The  solution  to  the  variational  problem  formulated  herein  predicts  the 
location  of  the  point  of  incipient  separation  xQ  on  the  airfoil  upper 
surface  and  the  magnitude  of  the  peak  velocity  (a  constant  from  the 
leading  edge  to  x0). 
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6.1  Small  Disturbance  Transonic  Assumption 

In  formulating  a  variational  problem  for  obtaining  the  maximum  lift 
achievable  by  a  transonic  airfoil  section,  use  has  been  made  of  the 
transonic  small  disturbance  approximation 

i(x,y)  *  U IJx  +  l2/3  *<x,y)  +  0(d4/3)] 

p(x,y)  *  p„[l  +  52/3  P(x,y)  +  0(64/3)l 

p (x,y)  *  pjl  +  <S2^3  o(x,y)  ♦  0(  54^3) ] 

-  /*5  * 

where  y  *  <S y  is  a  compressed  lateral  scale,  <5  is  the  thickness  ratio 
and  <)>,  P  and  a  are,  respectively,  the  perturbation  potential,  pressure 
and  density.  See  Figure  6. 

The  relationship  between  perturbation  pressure  and  velocity  is  given  by 

P  a  -  r  *x 

where  y  =  Cp/Cy  is  the  ratio  of  specific  heats. 

For  a  unit  chord  airfoil,  the  lift  coefficient  is  given  by 

cL  *  l/4  ».u  b 

where  l  »  f  Pw  «  -  f  Pupper  df 
Jo  Jo 

and  the  subscripts  "lower"  and  "upper"  refer  to  the  respective  sides  of 
the  airfoil . 
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Combining  the  previous  four  equations  gives 


where  =•  U^a,,.  is  the  freestream  Mach  number  and  a,,  =  yjy  p is  the 
corresponding  speed  of  sound.  Thus,  maximizing  lift  is  equivalent  to 
maximizing  M^  C.  .  To  the  leading  order  approximation,  the  transonic 
result  for  and  its  incompressible  analogue  for  CL  have  the  identi¬ 

cal  format.  Consequently,  the  same  ideas  developed  by  Liebeck  and  Smith 

[2]  to  obtain  the  optimum  (high  lift)  pressure  distribution  can  be  ap- 

2 

plied  to  the  transonic  case  as  well.  Therefore,  to  maximize  M^  CL,  we 
require  that 


on  the  lower  surface  be  an  absolute  minimum  while  on  the  upper  surface 
it  should  be  a  maximum.  If  ug/Ua#  could  be  made  to  vanish  along  the  en¬ 
tire  lower  surface  then  the  optimal  lower  surface  condition  would  be 
achieved.  However,  since  stagnation  can  only  occur  at  a  point  in  a  two- 
dimensional  flow,  a  more  realistic  requirement  is  that  I U0/UJ 1ower  be 
as  near  zero  as  possible.  In  general,  at  the  trailing  edge  u  =  ut  >  0 
Is  specified.  Since  flow  always  accelerates  from  a  stagnation  point, 
the  lower  surface  problem  seeks  the  airfoil  shape  that  permits  the  flow 
to  accelerate  away  from  the  stagnation  point  and  achieve  the  trailing 
edge  value  u^  in  such  a  way  that^ugd5  is  as  small  as  possible.  Since 
the  mathematical  limit  of  zero  area  under  the  curve  cannot  be  physically 


For  a  cusped  trailing  edge,  u^  tends  to  U„.  For  a  finite  T.E. 
angle,  the  Invlscid  u*  *  0.  The  ut  that  we  specify  corresponds  to 
the  tangential  flow  at  the  trailing  edge  external  to  the  boundary 
layer.  The  boundary  layer  tends  to  smooth  the  finite  T.E.  angle. 
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met,  the  problem  will  be  addressed  using  a  "cut  and  try"  method  investi¬ 
gating  the  various  practical  possibilities  for  achieving  the  constraint. 

The  upper  surface  problem  seeks  the  greatest  area  under  the  curve 
fQ^ ufid5  .  This  problem  Is  controlled  by  boundary  layer  separation  due 
to  adverse  pressure  gradient  in  the  recovery  region. 


6.2  Investigating  the  Quantity  to  be  Maximized 


We  now  focus  attention  on  the  upper  surface  problem  and  define  the  quan¬ 
tity  u/U.  that  Is  to  be  maximized. 


To  obtain  the  largest  value  of  H/U^  under  the  constraint  that  the  flow 
remains  fully  attached,  use  is  made  of  the  limiting  flow  on  the  verge  of 
separation.  As  has  been  developed  previously,  an  expression  for  the 
compressible  turbulent  boundary  layer  velocity  distribution  that  remains 
on  the  brink  of  separation  (in  terms  of  the  equivalent  flat  plate  length 
scale  s)  Is  given  there  by  [c.f.  equations  (5)  and  (6)]. 


ue(s) 


|  J£[l  -(b  (S/S0)1/6-  1)1/3] 

£  (w  *  °„)'1/4 


1/2 


•*  so  i  5  i  sc 


5  SC  1  5  -St 


(22) 


where 


a5  3  °3  so 


uQ  Is  the  peak  velocity  attained  at  s  =»  sQ,  st  denotes  the  equivalent 

flat  plate  trailing  edge  length,  s  denotes  the  location  along  an  equlv- 

★  C 

alent  flat  plate  when  Cp  *  4/7,  and  B,  o^,  are  specified  functions  of 
the  peak  Mach  number  ^  and  equivalent  flat  plate  Reynolds  ntrober  R  . 

Figure  7  Illustrates  the  linear  relationship  between  the  airfoil  arc 
length  measured  from  the  nose  x  and  the  equivalent  flat  plate  length 
scale  s  along  with  definitions  of  the  constants  k  and  K  (the  two  key 
parameters  of  the  variational  procedure  to  follow).  That  is, 

5'so  *  *-*„  : 

k  i  *0/s0  (  (23) 

*  5  VSo 

The  Integral  appearing  In  equation  (21)  can  be  expressed  as  a  contribu¬ 

tion  from  the  accelerating  region  and  a  contribution  from  the  pressure 
recovery  region.  Breaking  up  the  integral  in  this  way  and  using  the 
linear  relationship  between  s  and  x  of  equation  (23)  in  the  pressure 
recovery  region  only,  gives 


Acceleration  Region  Pressure  Recovery 


*  Note  we  have  arranged  that  In  the  acceleration  region  the  integrand  is 
a  function  of  the  physical  arc-length  [i.e.,  u  =  u(0]  and  in  the 
pressure  recovery  region,  the  Integrand  is  a  function  of  the  equiva¬ 
lent  flat  plate  length  s. 
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Si  nee  the  velocity  distribution  in  the  pressure  recovery  region  is  spec¬ 
ified  by  equation  (22),  the  second  integral  in  the  above  expression  can 
be  written  as 


IS(K,RS  ,M  ) 


/•< 

A 


u  (s) 

D —  d(s/so) 
00 


(24) 


and  therefore 


u 

v: 


IJK.R  ,M  ) 
o 


d(C/s0) 


(k  +  K  -  1) 


(25) 


The  Reynolds  number  and  Mach  nianber  functional  dependence  enters  into 
the  Is  expression  by  way  of  equation  (22).  It  can  be  shown  that  the 
peak  Mach  number  M0  depends  only  on  uQ  for  specified  freestream  condi¬ 
tions  since 

u  2  a2  U2  a2 

o  .  o  -  .  « 

~7T  ITTT 3  ~T~  TTTT 


Hence, 


U.  a 
00  00 

,T"  +  T7TT 


-V) 


(  Y-l) 


(26) 


It  can  also  be  shown  that  the  equivalent  flat  plate  Reynolds  number 

R  depends  functionally  only  on  u0  and  s  provided  the  freestream 
so  0 

conditions  are  known.  Suppose  the  freestream  unit  length  Reynolds  num¬ 
ber  is  R.  =  Then 


(27) 
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where 


v-  u-  p0 

V  *  U  P„ 

0  0  * 

Using  the  power  law  relationship  between  u  and  T  (see  Reference  [10], 
NACA  Report  1135) 


and  the  isentropic  relationship 


R. 


fc) 


1/(Y-1) 


Y-l 

~T~ 

~T~ 


1 

TTTF 


in  equation  (27)  gives 


-  K 


1 1  *  ^  «.2 


(7-3y) 

TITTT 


(28) 


which  displays  the  u0  and  s0  dependence.  Furthermore,  uQ  and  s0  depend 
on  the  parameters  k  and  K  since  from  equation  (23) 


sQ  »  l/(K+k-l)  j  (29) 

xQ  *  k/(K+k-l)  ) 


and  from  equation  (22) 


u0/U.  -  f(k,K;  ut/Uj 


(30) 


In  equation  (30),  the  trailing  edge  velocity  ratio  ut/Us>>  is  assumed 
given  for  our  purposes. 


As  a  consequence  of  equations  (26)  through  (30)  we  have  determined  that 
the  peak  Mach  number  and  the  flat  plate  Reynolds  number  are  functions 
only  of  k,  and  K  provided  the  freestream  conditions  and  trailing  edge 
velocity  ratio  u^/U,,  are  specified  in  advance.  Therefore,  In  equation 
(25),  I  (K,R  ,M  )  =  J(k ,K) . 

5  •>«  u 


6.3  Two  Limiting  Cases  of  the  Equivalent  Flat  Plate  Length 

The  two  special  cases  to  be  considered  for  this  investigation  are  speci¬ 
fied  by  xtr*x0  (full  laminar  run  to  the  separation  point)  and  xtr=0 
(full  turbulent  boundary  layer).  For  the  laminar  run  case  utr=u0, 
vtr’vo  and  e<luation  (20)  reduces  to 

3/5 

x  ;  laminar  run  (31) 

o 

For  a  full  turbulent  boundary  layer  (x^r=  0  ;  9trs  0),  equation  (30)  re¬ 
duces  to 


s  * 


1.71 


3-94f1/U JO  ^3.94 


Re 


1/25 


[ft) '  /  © 


a(5/x0) 


24/25 

x 


o’ 


full  turbulent  (32) 


Since  k  s  x0/s0  .  equations  (31)  and  (32)  provide  definitions  of  k  for 
the  laminar  run  and  full  turbulent  cases,  respectively.  That  is 
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fY< 

w. 

-  ©V3/ 

1  kA2/3  (970.7)  =  0 

(33) 

JQ 

A 

k 

»\3-94  /uo\3‘94 

:)  d(5/so)  •  (vz) 

Rx1/2Vl.749  kt1/24  -  0 

0 

(34) 

Jo 

u  X 

where  Re  = 
xo 

-2-2-  and  kA,kt  denote  k 

0 

for  the  laminar  run  and 

turbul ent 

cases,  respectively. 


6.4  Variational  Formulation  and  Solution 

For  both  the  laminar  run  and  turbulent  cases,  the  variational  problem 
can  be  formally  written  by 

rk 

f[Y(n)]dn  ♦  J(k,K) 

(maximize)  HCY(n) ;k,K]  *  — -  (35) 

G20c,K1 

under  the  constraint 
>k 

G(k,K)  +  g[Y(n)]dn  -  0  (36) 

^o 

That  is,  it  Is  desired  to  find  the  function  Y(n)  and  the  values  of  k  and 
K  which  maximize  the  functional  H[Y(n);k,K]  subject  to  a  side  condition, 
equation  (36).  Equations  (33)  and  (34)  are  the  side  conditions  for  the 
two  cases  of  Interest.  Compare  also  equations  (25)  and  (35). 


Introducing  the  Lagrange  multiplier  A,  equation  (35)  and  (36)  may  be 
combined  in  the  following  way: 


H[Y(n);k,K,A] 


I[Y(n)  ,k;  A]  ♦  Gjtk.K.A) 


(37) 
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where 


I[Y(n)k;X] 


G^k.K.X)  ■  J(k,K)  +  X  G(k,K) 


(38) 


The  first  variation  of  H  is 

rij  3H  »  .  3H  .  »  .  3H  *  /* 

6H  3  IT  51  +  SG^  Gl  +  ^  G2 

Each  of  the  variations  61,  6G^,  SGg  can  be  expressed  in  terms  of  the 
variations  6Y,  6k  and  6K.  Letting  6H  *  0  gives  the  following  Euler 
equations 


f  +  X  g 


5Y  *  0 


3H 


ij  (f  +  Xg) 


rpk 


3H  3G1  3H  3G2  .  - 

w +  nq  -sir  0 


3H  3G1  3H  3G2  3  n 

^  nr +  -ssj  nr  0 


6k  *  0 


5K  *  0 


(39) 

(40) 

(41) 


Consider  first  the  Euler  equation  (39)  to  obtain  the  shape  of  the  velo¬ 
city  profile  in  the  region  0  <_  x/sQ  <_  k.  From  equation  (25) 


*  a  VU- 

f[Y]  =  Y 


(42) 


For  the  laminar  run  case,  equation  (33)  gives 
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g[Y]  =  r 


(43) 


Therefore,  substituting  (42)  and  (43)  into  (39)  gives 
.  /  u^ 


1  +  5  X 


;  laminar  case 


(44) 


The  6Y  Euler  equation  gives  for  the  turbulent  case  (using  equation  (34) 
from  which  g[Y]  »  YB*94): 


1  +  3 


I  ue\ 
.94  \\rf-) 


u_\2.94 


;  turbulent  case 


(45) 


Both  equations  (44)  and  (45)  imply  that  ue/U.  is  constant  in  the  range 


0  i  x/s0  i  k*  But  Ue/U«  “  uo/U-  at  xo/so  3  k  »  so 


VU-  s  VU- 


0  <_  x/sQ  <_  1c 


(46) 


The  corresponding  Lagrange  multipliers  can  be  obtained  from  equation 
(44)  to  (46). 


Solving  equations  (40)  and  (41)  provide  the  value  of  k  and  of  K  which 
maximize  TT/U^.  Figure  8  illustrates  the  physical  situation  where  the 
particular  k  and  K  so  obtained  fix  the  peak  velocity  u0  and  its  location 
xQ  on  the  airfoil.  It  is  difficult  to  obtain  k,K  analytically  due  to 
algebraic  complexities.  Consequently,  a  numerical  procedure  was  used. 
Figure  9  and  10  display  the  variation  of  maximum  theoretical  lift 
coefficient  (for  both  the  laminar  run  and  fully  turbulent  cases, 
respectively)  with  for  a  range  of  u^/U,,.  Figures  11  and  12  present 
similar  results  for  maximum  theoretical  C,  .  It  is  seen  that  a 
substantial  advantage  exists  for  the  laminar  run  case.  Furthermore,  the 
sensitivity  of  the  results  to  the  parameter  ut/U„  is  quite  apparent. 
Structural  considerations  at  the  trailing  edge  of  an  airfoil  apparently 
1  imit  UJ./U,,  to  about  0.9  since  our  perusal  of  airfoil  data  has  not 
displayed  values  any  greater. 


u 


Family  of  Ncnseoarati nc  “1st  Fooftoo  V 
Distributions  for  i  Given  u„/’j  zr.c  M_ 


ricure  8 
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6.5  Modifications  To  The  Optimized  Velocity  Distribution 

The  analysis  of  the  previous  sections  has  Indicated  that  CL  will  be  max¬ 
imized  by  using  a  velocity  distribution  of  the  form: 

ue 

— s  o  over  the  entire  lower  surface 

^ao 

—  given  by  a  flat  roof-top  plus  compressible 
u»  Stratford  distribution  on  the  upper  surface. 

Unfortunately,  this  distribution  will  not  yield  a  physically  realizable 
airfoil  shape  due  to  the  discontinuities  present  at  the  nose  and  the 
peak  velocity  point  and  the  fact  that  true  stagnation  can  not  occur  over 
the  entire  lower  surface.  Therefore,  It  seems  necessary  to  modify  the 
velocity  distribution  around  the  nose  region,  over  the  entire  lower  sur¬ 
face,  and  at  the  sudden  decelerating  region  located  in  the  supersonic 
flow  on  the  upper  surface  {x=x0)  to  avoid  a  possible  shock  there.  The 
proposed  velocity  distribution  modifications  are  indicated  by  the  broken 
lines  in  Figure  13.  Systematic  variation  remains  yet  to  determine  how 
close  to  the  optimum  pressure  distribution  can  be  obtained. 
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7.  SHOCK  CONSIDERATIONS 

Our  design  procedure  relies  on  relating  the  developing  boundary  layer 
along  an  airfoil  surface  with  that  on  an  "equivalent"  flat  plate  (the 
equivalence  determined  via  the  momentum  thickness).  Consequently,  it  is 
important  for  our  design  procedure  to  have  the  capability  of  assessing 
the  shock -boundary  layer  interaction  through  its  effect  on  the  momentum 
thickness.  Appendix  C  presents  a  derivation  of  the  appropriate  expres¬ 
sions  necessary  to  assess  the  shock-boundary  layer  Interaction  for 
Incorporating  the  shock  effect  into  our  design  procedure. 

The  momentum  thickness  after  the  shock  9$  was  found  to  be  related  to  the 
momentum  thickness  before  the  shock  by 

9  =  F(M  ;k)  9  (47) 

S  ^  0  0 

where 

F(Mq ;k J  =  1 

A 

g  *  0.713  is  a  compressibility  factor  for  0.8  <  M  <  1.4 


v  Q.:.llk+Ulk4£l  U+4MJU-M'2) 
L  (7-1)  k  g  J  0  0 
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where  is  the  Mach  number  just  downstream  of  the  shock.  M^  is  related 
to  M0  by  the  normal  shock  jump  relations. 

A 

For  typical  values  of  k  (-6)  and  MQ  (-1.2),  the  change  in  the  momentum 
thickness  across  the  shock  (relative  to  the  momentum  thickness  before 
the  shock)  is  0(1).  Equation  (18)  provides  an  expression  for  aQ  for  the 
full  laminar  run  case  and  the  full  turbulent  case. 


To  obtain  the  equivalent  shock  free  flat  plate  length  sso  (comparable  to 
s0  in  the  case  of  no  shocks),  the  following  expression  is  used: 


The  above  equation  provides  an  implicit  relationship  between  the  equiv¬ 
alent  flat  plate  length  and  the  shock  location  x0  on  the  airfoil. 


Now,  if  Cp 


★ 


(s) 


is  defined  by 


(50) 


for  the  airfoil  having  a  shock  at  xQ  and  a  velocity  immediately  down¬ 
stream  of  the  shock,  then  the  separation  criterion  is  given  precisely  by 
equation  (2)  but  with  n  determined  by  equation  (48)  and  sQ  replaced  by 


In  this  case,  Cp  for  Incipient  separation  is  given  by 


*  n  s  v1/16  1 

vsi  ■  to  -iJ  ■■ 

Cp(s)  »  1  -  'o3s*V-1/2; 


siSsc 


Si$SC 


(51) 
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From  this  last  equation  2  <  n  <  k. 


7.1  Variational  Problem  Including  the  Shock 

The  variational  problem  including  a  normal  shock  follows  along  the  same 
lines  as  in  the  shockless  case  of  Section  6.  The  parameters  k  and  K 
(see  Figure  14)  are  defined  for  the  shock  case  by 


k  *  Vsso 


1-x. 


'so 


+  1 


and  the  integral  Is  (see  equation  (24))  Is  given  by 

!s(Ki  k  *  vMi) H  f+  ^ 

•  o  Jj  +  •  so ' 
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where 


U1 

c 


U1 
I r 


if 


K 


<  K 


The  location  of  the  shock  sso  is  obtained  directly  from  the  expression 
9s  3  POV^o  9'iven  by  equations  (47)  and  (C-1S).  Using  the  defini¬ 
tion  of  ksx  /s_  and  the  above  mentioned  equations  leads  to: 
o  so 


5/24  -1/4 

£)3'9V*)*  °-5725  k’I/24  Hn°*rS'\*k)  fe)  ,521 


By  using  the  conditions  along  a  streamline  (i.e.,  enthalpy  remains 
invariant,  jump  conditions  across  the  shock  and  isentropic  flow  else¬ 
where)  It  can  be  shown  that  the  right  hand  side  of  equation  (52)  depends 
functionally  upon  k,  K  and  parametrically  upon  M^.v^  and  the  speed  ratic 
at  the  trailing  edge  u^/U^,  Furthermore,  the  variational  formulation 
for  the  shocked  flow  case  can  be  shown  to  reduce  to  exactly  the  form  of 
the  shockless  flow  case  considered  previously  with  the  result  that 

V  uo  ;  (°ixixo=k  sso^  (53) 

Consequently,  the  Integral  expression  on  the  left  hand  side  of  equation 
(52)  Is  equal  to  k.  Therefore,  equation  (52)  provides  an  implicit  rela¬ 
tionship  between  k  and  K.  That  is,  equation  (52)  is  of  the  form  k=k(K) 
and  xQ,  sSQ,  uQ  andf  ue(x)/U„  dx  depend  functionally  only  on  K  and 

parametrically  on  and  u^/U,,.  To  determine  the  maximum  lift  a 

numerical  procedure  was  developed. 


Numerical  Results  for  the  Maximum  Lift  Problem 

In  the  shocked  flow  problem  the  order  of  n  (the  exponent  of  the  power 
law  of  the  velocity  profile)  is  determined  from 

M  2 

(  (k+3  j  (/)  -  3  ;  Mj  <  1 

"  a  \  0 

^  k  ;  _>  1  (shock  free) 

It  was  determined  that  the  shock  solutions  matched  the  shock  free  solu¬ 
tions  only  when  uQ  *  aQ  (Mq*  l;K=Kson.|c  j.  See  figure  15.  In  order  to 
calculate  the  shock  solutions,  an  iteration  method  was  employed  using 
the  exact  shock  free  sonic  solution  as  an  initial  guess.  Some  results 
are  provided  in  Figures  16  -  18  for  the  trailing  edge  speed  ratio 
ut/U„  »  0.9  and  for  a  range  of  free  stream  Mach  numbers  M,,,  *  0.8,  0.9 
and  1.0.  As  M*  Increases,  the  length  of  the  roof-top  pressure  distribu¬ 
tion  Is  seen  to  increase.  Furthermore,  the  Stratford  type  velocity 
profile  exhibits  less  streamwise  variation  as  the  strength  of  the  shock 
Increases.  This  is  expected  on  physical  grounds  since  the  effect  of  an 
Increase  In  shock  strength  is  to  bring  the  boundary  layer  velocity 
profile  closer  to  the  condition  of  separation.  Figure  19  shows  for  the 
trailing  edge  speed  ratio  u^/U.  *  0.9  how  the  maximum  lift  coefficient 
was  found  to  vary  with  free  stream  Mach  number  in  the  transonic  speed 
regime  for  a  fully  turbulent  boundary  layer  for  the  cases  of  shocked  and 
shock  free  flow.  The  minimum  percentage  difference  between  the  case 
with  shocks  and  without  shocks  was  found  to  be  approximately  11%;  with 
the  occurrence  of  shocks  diminishing  the  maximun  lift  coefficient. 


Figure  18.  A  Family  of  Shock  Solutions  for  Given  Values  of  and  u./U 
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Figure  19.  Maximum  CL  With  and  Without  Shock 
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I nner  Field  Development 

The  Inner  field  region  corresponds  to  the  flow  within  the  turbulent 
boundary  layer  near  the  wall  surface.  Within  this  region,  we  assume 
fully  turbulent  flow  and  use  the  Prandtl  mixing  length  hypothesis  to 
represent  the  shear  stress  with  A=<y  where  <  is  von  Karman's  constant 
and  y  is  the  local  measure  of  transverse  distance  along  the  surface. 
The  shear  stress  is  given  by  the  expression 

T  *  p*V  (17)  (A~U 

By  integrating  the  boundary  layer  equations  from  the  wall  outward  to 
some  point  y,  which  is  also  near  the  wall,  we  obtain 

T  =  Tw  +  y  TJs  ;  P  =  P<^  (A-2) 

where  the  external  pressure  gradient  has  been  assumed  to  be  impressed 
upon  the  boundary  layer. 

Equating  t's  In  expressions  (A-l)  and  (A-2)  gives  the  flow  speed  gra¬ 
dient  at  a  point  y*  where  the  two  expressions  for  t  are  equal.  That  is. 


Assuming  separation  occurs  when  xw  vanishes,  expression  (A-3)  can  be 
viewed  as  a  partial  differential  equation  for  the  flow  speed  in  a  region 
of  sufficiently  small  y  .  That  Is,  setting  xw  =  0  in  the  previous  equa¬ 
tion,  gives 

1/2 

(y  |y  3  A  y'1/2  when  xw  =  0  (A-4) 

where 


c  =  u/uQ  ,  a  speed  ratio  scaled  with  the  reference  speed  uQ. 

(15  9  \  1  /  2 

p"  k’^u”  -»§■  ,  a  pressure  gradient  coefficient 

’  (varying  with  location  s) 


Before  equation  (A-4)  can  be  Integrated,  we  must  determine  how  the  den¬ 
sity  ratio  p/pw  depends  on  the  streamwlse  speed.  For  a  perfect  gas,  the 
density  ratio  Is  related  to  the  temperature  ratio  by 


From  Crocco's  relationship  (see  Reference  [9],  page  144),  we  have  that 


(1  +  -  A2?2] 


where 


A 


2 


[( y-1)/2]  u2 

T7T —  Me 
w  e 


1  +  [( y-1)/2]  M2 

bi  - 


(A -6) 


(A-7) 


Therefore,  combining  equations  (A-4)  -  (A-7)  gives 

-1/2 


H 

3y 


Ch-Q1«;-a2c2  )1/2 


A  y 


when  t  »  0 
w 


( A— 8 ) 


This  equation  can  be  integrated  as  an  ordinary  differential  equation 
since  s  enters  only  as  a  parameter  through  the  terms  A,  3^  and  A.  Inte¬ 
grating  from  the  wall  outward  the  solving  for  c  as  a  function  of  y 
yields 

C  *  5  sin  (2AA^y-c)  +  n 


( A-9) 
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where  . 

i  5  (4A2  +  B2)  /2A2  1 

n  5  Bj/ZA2  >  ( A— 10) 

e  s  sirf^/UA2  +  B2]1^  J 

To  approximate  the  flow  speed  ratio  near  the  wall  at  separation,  we 
expand  the  sine  function  appearing  in  equation  (A-9)  using  the  trigono¬ 
metric  relation  for  the  sine  of  the  difference  of  two  angles  together 

with  the  series  expansions  for  the  sine  and  cosine  of  small  angles. 

Neglecting  higher  order  terms  in  the  expansion,  we  obtain 

(  y  «  1 

C  •  2A  y/y  [l  +  £  BjVJ*  0(y>]  ;  {  (A-U) 

(  Tw  =  0 

Notice  that  the  flow  speed  tends  to  zero  in  a  square  root  fashion  at 

separation  as  the  wall  is  approached.  The  leading  order  term  is  similar 
to  the  incompressible  result  (suitably  modified  by  the  compressible 
pressure  gradient  term  In  A)  and  the  next  term  contains  a  heat  transfer 
contribution  and  a  Mach  number  effect  (see  equation  (A-7)l.  The  stream- 
wise  speed  gradient  at  separation  is  obtained  by  differentiation  of 
equation  (A-ll).  Hence, 

r  (  y  «  1 

•  A  y'1/2  [l  t  A8jy1/2  *  0ly>]  ;  J  iA-12) 

(  -  0 

and  the  slope  of  the  flow  speed  is  infinite  at  separation,  tending  there 
as  the  Inverse  square  root  of  the  transverse  distance. 
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To  complete  the  near  field  specification,  we  require  the  stream  function 
in  the  boundary  layer.  A  suitable  definition  of  the  stream  function  for 
compressible  flow  is 


But  for  y  <<  1  (which  implies  that  c  <<  1  from  equation  ( A- ID) 


p 


(1  +  B1c  -  A2c2} 


1  -  B1«  +  0(;2) 


(A-14) 


Substituting  equation  (A-ll)  Into  equation  (A-14)  and  the  resulting 
expression  into  equation  (A-13)  gives,  upon  integration 


♦ 


4 

7 


A  y 


3/zri ♦ 1 


L 


AB.y 


1/2 


j  y  «  l 

l  \ s  ° 


(A-15) 


Later  in  the  development  we  will  equate  certain  inner  and  outer  physical 
quantities  in  an  overlapping  region  to  obtain  the  separation  criter 
ion.  The  technique  is  not  an  "asymptotic  matching"  in  the  strict  sense, 
but  rather  falls  into  the  category  of  a  "patching"  approach.  The  proce¬ 
dure  Is  not  unique.  Nevertheless,  we  follow  Stratford  in  this  develop¬ 
ment  in  view  of  the  success  he  obtained.  In  so  doing,  we  require  the 

limiting  forms  (for  y  <<  1  and  tw  =  0)  of  the  quantities  and 

These  are  readily  obtained  from  equations  (A-ll),  (A-12)  and  (A-15). 
They  are 
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3  (  y  «  i 

*(|f)  3  7  uo  ^L1  +  I  ABlV*  +  0(y)]  ;  {  (A-16) 


U’° 

The  physical  significance  of  these  quantities  is  unclear.  However,  it 
is  perhaps  interesting  that  equations  (A-16)  and  (A-17)  have  the  expres¬ 
sions  In  brackets  in  common  -  at  least  to  0(y}.  The  ratios  were  origin¬ 
ally  chosen  (most  likely)  to  clear  the  y  dependence  from  the  leading 
order  terms. 


Outer  Field  Development 

The  outer  field  is  considered  to  be  the  region  in  the  turbulent  boundary 
layer.,  well  removed  from  the  surface  of  the  body.  If  a  pressure  rise 
did  not  exist  in  this  region,  the  velocity  profile  would  be  similar  to 
that  along  a  turbulent  flat  plate.  The  total  pressure  loss  can  be  cal¬ 
culated  for  such  a  case.  For  the  real  case  of  a  pressure  rise  we  make 
the  assumption  that  the  loss  of  total  pressure  along  a  stream  line  is 
Independent  of  the  pressure  rise.  This  assumption  is  supported  by 
experimental  results.  (See  Reference  [3],  Grabowski  st  al.).  If  the 
pressure  rise  is  severe,  the  assumption  should  be  even  more  reliable 
since  it  will  produce  the  effect  of  an  instantaneous  change  of  flow  con¬ 
ditions  over  a  short  spatial  duration  and  for  such  a  change  the  effects 
of  dissipation  are  negligible.  As  a  consequence,  we  assume  that  the 
pressure  loss  along  the  stream  line  in  the  actual  case  is  the  same  as  on 
a  corresponding  stream  line  in  the  turbulent  flat  plate  case  where  the 
pressure  is  constant. 


In  the  outer  region,  we  assume  the  Bernoulli  equation,  applied  along  a 
stream  line,  holds  true  in  the  form: 


( A— 18) 

where  «|>^  denotes  a  limiting  stream  line  on  the  inner  periphery  of  the 
outer  region.  is  unknown  and  can  be  found  as  part  of  the  solution.) 


u2 

T 


I 


d£ 

P 


f(  *) 


along  \i>  > 


In  this  outer  region  we  assume  the  flow  is  well  represented  by  an  adia¬ 
batic,  inviscid  and  nonconducting  fluid  flow.  Such  a  flow  is  isentropic 
and  is  described  by  the  isentropic  relation 


ik)  ■  iy 


vy/(y-1) 


( A-19) 


Substituting  Equation  (A-19)  into  the  Bernoulli  equation  yields 


y  po  r-/_  i(y-1)/y_  r/.,A 
rFTTlT^p/poJ  *  fU) 


'P  >_  ^ 


( A-20) 


Mow  assume  that  far  from  the  wall  a  comparison  velocity  profile 
uci  uc( s , 4»)  exits  (subscript  c  denotes  the  comparison  flow)  and  is  a 
reasonably  accurate  representation  of  the  flow  field  in  the  outer  bound¬ 
ary  layer.  (Notice  that  the  y  dependence  enters  implicitly  by  way  of 
1>.)  At  some  point  s  =»  s0  corresponding  to  the  beginning  of  the  pressure 
recovery  (the  reference  point),  we  have  that 


u(So.'f)  -  uc(s0,*J 


p  -  pr 


(A-21) 


The  relations  appearing  In  equation  (A-21)  are  used  to  evaluate  the 

Bernoulli  'constant'  of  equation  (A-20).  Dividing  equation  (A-20) 
2 

by  Ug/2  results  in 
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C2  *  il  -  C*(s) 
c  p 


where 


c;<s>  ‘  ty^tt  - 


TVo 


(!U 


Yii 

Y 


-  1 


ujs,^) 

Cc  5  -V - 


C  5 


u(S,! 


( A-22) 


Taking  the  partial  derivative  with  respect  to  ♦  of  the  Bernoulli  inte¬ 
gral  (A-22)  gives 


-  3C  .  3?C 

c  7?  '  Cc  W~ 


( A—  23) 


From  the  definition  of  the  stream  function  given  by  equation  (A-13)  we 
have  that 


%  ■  “o (k> 


[%\  ■  '« 


Therefore, 


IT  *  (ly)c  ly13  uo  (y  ?c 


( A— 24 ) 


3?, 


Comparing  Equations  (A-23)  and  (A-24)  and  assuming  (p/p  )  -  (p/p  ),  we 
obtain  c 


>c  ,  3!c 

?  If 

Using  Crocco's  relationship  for  the  comparison  Flow  we  obtain 


( A-25) 
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iy 


=  (l  +  b1cc-a2s|) 


4  (i  ♦ 


'»  t1  *  bi«c-*2!c> 


dy 


3C  .  3CC 

W 

and  from  equation  (A-22) 


C2  3  4  -  C*(s) 


(A-26) 


(A-27) 


( A-28) 


( A— 29 ) 


Equations  (A-29),  (A-28)  and  (A-27),  representative  of  the  outer  field, 
are  comparable  to  the  inner  field  equations  (A-ll),  (A-12)  and  (A-15), 
respectively. 

We  now  assume  the  following  power  law  expression*  for  the  comparison 
flow: 

Cc  *  Aym  ;  A  =  (A-30) 


where  the  power  m  of  y  and  the  coefficient  A  (expressed  in  terms  of  a 
constant  k^  determined  from  experimental  data  and  the  momentum  thickness 
9)  are  given  by  turbulent  flat  plate  boundary  layer  experiment.  Experi¬ 
mental  evidence  exists  showing  m  *  1/7  fits  data  well  even  for  a  Mach 
number  as  great  as  2.4. 

Equation  (A-30)  is  to  be  substituted  into  equations  for  c  and  3c/ 3y 
(previously  given).  The  immediate  objective  is  to  develop  the  far  field 
patching  conditions.  In  so  doing,  we  will  push  the  far  field  results  to 


*  Other  velocity  profiles  can  be  chosen,  for  example,  Cole's  law  of  the 
wall /I aw  of  the  wake.  However,  the  assumptions  inherent  in  this 
Stratford  type  development  do  not  warrant  such  sophistication. 
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the  limit  of  small  y.  Patching  will  then  provide  y  at  the  join  (yayj) 
and  a  relationship  of  the  actual  velocity  profile  to  the  comparison  pro¬ 
file.  Hence,  for  small  y 


a?  .  3cc 
ly  Jy~ 


\  m  y" 


y  «  1 


Jry 

|  cc(l-8l!c)dy  ■  V  ferr  [l-hx  lS?TT  ym];  * «  1 

A 


( A—  31) 


The  patching  conditions  are  readily  obtained  from  these  expressions. 
They  are 


u  x4  A-rr  y41"'2 
o  (m-1)  3 


m+1 

IShT 


V  *"] 


(A-32) 


Notice  that  the  expresions  in  brackets  above  are  not  identical,  as  they 
were  found  to  be  for  the  similar  conditions  obtained  in  the  inner 
field.  This  disparity  arises  because  the  approximation  for  the  stream 
function,  flow  speed  and  transverse  flow  speed  gradient  (using  the  com¬ 
parison  profile)  become  weak  for  y  «  1.  However,  to  leading  order  in 
y,  this  disparity  does  not  upset  the  approach  we  are  using.  It  would, 
however,  affect  matters  if  a  higher  order  development  were  pursued.  In 
that  case,  a  better  representation  of  the  velocity  profile  in  the  outer 
boundary  layer  for  small  y  would  be  warranted. 
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Furthermore,  it  should  be  noted  that  the  relationship  between  c,  c 
*  c 

and  Cp  (equation  ( A-29 ) )  was  not  used  in  obtaining  equation  (A-33). 
After  equating  the  near  field  and  far  field  results,  equation  (A-29) 
provides  the  separation  criterion. 


Patching 


To  find  y  at  the  patching  locating  (y»yj)  we  equate  the  coefficients  of 
equations  (A-16)  and  (A-32)  and  obtain 


4 

T 


m3 

Tm+TT 


4m-2 


( A— 34 ) 


Solving  this  expression  for  y^/e  (see  equation  (A-30) )  gives 

4 


(y9) 


2 -4m  3  m4  k'l 


Tii+n 


e2  a4 


( A-35) 


To  find  the  relationship  between  the  actual  velocity  profile  and  the 
comparison  profile  at  the  patching  location,  we  equate  coefficients  of 
equations  ( A- 1 7 )  and  (A-33)  and  obtain* 


( A— 36 ) 


*  Notice  that  3  ^y^Cl-Sm/tm+l)]  - 

C*  <  (l-2m)/( 1+m)  ^  (n-2)/(n+l)  ;  m  =l/n 

P  - 
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Separation  Criterion 


Equations  (A-29),  (A-35)  and  { A-36 )  contain  the  essential  results  of  the 

analysis  thus  far.  Since  we  seek  a  separation  criterion  in  terms 
★ 

of  C  ,  it  is  necessary  to  express  dp/ds  appearing  in  the  A  expression 
P  ★ 

(see  A  definition  occurring  after  equation  (A-4))  in  terms  of  Cp. 


From  equation  (A-22)  and  the  isentropic  gas  relationship  (A-19) 


(A-37) 


It  should  be  noticed  that  the  density  pg  appears  in  the  pressure  gra¬ 
dient  expression  rather  than  the  reference  density  p0  or  pw. 


,4a2  u 
* 


Consequently,  the  term  in  Equation  (A-35)  becomes 

2 

dC.' 


aV 


and  equation  (A-35)  becomes 


(>y9) 


2 -4m 


{£)  fe)  (• S) 2 


( A— 38 ) 


( A-39) 


The  density  ratio  can  be  expressed  as 


1 


( l+Bj-A2)  w 


( A-40) 


With  the  additional  assumption  that  the  wall  is  an  adiabatic  wall 

(vTJ  °“"- 


( A— 41 ) 


and  since  the  local  Invlscid  Mach  number  varies  along  the  airfoil  sur¬ 
face  the  quantity  (pw/pJ  is  a  function  of  s. 

w  e  a 

The  recovery  factor  r  is  approximately  equal  to  0.9  (see  Reference 
[9]).  Equation  (A-39)  can  be  written  as 


1/4  „  1/2 


*  -1/2 


(£)  (•?) 


( A-42 ) 


Combining  Equations  (A-29),  (A-30)  and  (A-36)  and  solving 

for  (y./a)2111  gives 


;y  /9l2n»  =  P.  -  .  •  r3n» 

0  2  -  *  1  '  ^ 
2  Ml1  WITj 


{ A-43 ) 


Eliminating  (y ^/Q ]  from  equations  (A-42)  and  (A-43)  gives  the  separation 
criterion.  Letting  m  =  1/n  recovers  Stratford's  form  of  the  criter- 
i on .  That  i s , 

^(n-2)/  do*/72  "r  -11/4  /n.2^"'2Vw\1/2 

[s  (e  3/)  -  *1  ~r~  ST)  $  tA-44' 

Ln^(n+l)J  e  a 


where  0  £  C  <  (n-2)/(n+l)  is  required  (see  footnote  on  page  66). 
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Compressible  Momentun  Thickness 


Stratford  expresses  0  as  a  function  of  local  Reynolds  number  Rs  and  s 
for  the  incompressible  case.  One  way  to  proceed  is  to  look  at  the  local 
skin  friction  coefficient  c^  vs.  Me  curves  as  in  Clutter's  report 
[11].  These  are  based  on  the  Van  Driest  II  formula  and  are  considered 
state  of  the  technology. 


One  starts  with  the  momentum  integral  equation  for  a  flat-plate  boundary 
layer 


d8 

ls  = 


where 


T  * 


W 


( A— 45 ) 


If  a  power  law  form  for  cf  is  assumed 


then  equation  (A-45)  can  be  integrated  to  give  0  with  the  desired  func¬ 
tional  dependence.  For  Mach  numbers  in  the  range  0  to  2  and  Reynolds 
numbers  in  the  range  IQ6  -  107  the  variation  of  the  constants  a,b  were 
found  [4]  to  be  about  ten  percent.  Considering  the  approximate  nature 
of  the  Stratford  method  to  begin  with,  we  take  a,b  to  be  values  found 
when  Me  is  near  unity.  That  is, 


a  5  0.036 
b  5  -  1/6 


Then 


0  -  0.022  s  RJ1/6 
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The  separation  criterion  can  be  written,  therefore,  as 


This  is  the  form  of  the  separation  criterion  reported  in  the  introduc¬ 
tion  to  this  section.  The  density  ratio  on  the  right  hand  side  is  given 
by  equation  (A-41). 

From  Reference  [8]  (Schl ichting,  pg.  674,  giving  measurements  on  a  flat 
plate  at  zero  incidence  to  the  flow)  we  find  that  the  l/7th  power  law 
fitting  experimental  data  quite  well  even  when  Mach  number  equals  2.4. 
In  this  case,  the  empirical  constant  k^  is 

kj  *  0.683 

which  will  be  assumed  to  remain  valid  throughout  this  development  for 
other  values  of  n,  as  well. 

Stratford  [1]  expresses  his  separation  criterion  in  terms  of  the  pres¬ 
sure  coefficient  referenced  to  the  dynamic  head  at  the  peak  velocity 
point.  That  is. 


★  _ 

C„  i s  rel ated  to  C„  by 
P  P 


(A-47) 
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( A-48 ) 


From  this  last  expression  it  can  be  seen  that  C  *  Z  as  M  ♦  0  . 

p  p  o 


Another  way  of  expressing  Cp  is  in  terms  of  the  flow  speed  ratio 
ue/uQ.  From  Reference  [12]  (page  55)  we  have  that 
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(A-49) 


Combining  this  last  result  with  the  last  of  equation  (A-48)  gives 


( A-50) 


Since  the  pressure  coefficient  is  determined  from  conditions  prevailing 

external  to  the  boundary  layer  (the  inviscid  solution),  the  previous 
★ 

expression  for  Cp  could  have  been  deduced  immediately  from  the  Bernoulli 
equation  (A-22).  In  that  expression,  ;  *  1  along  the  edge  of  the  bound¬ 
ary  layer. 
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APPEMDIX  B 


Derivation  of  the  Compressible  Stratford  Flows 
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Compressible  Stratford  Flows 


The  density  ratio  on  the  right  hand  side  of  equation  (A-46)  depends  on 
the  boundary  layer  edge  Mach  number  Me,  as  given  by  equation  (A-41). 
From  the  Bernoulli  equation,  the  edge  Mach  number  is  related  to  Cp  and 


the  peak  Mach  number  M0  by 

/(1+IT-m2oC'p) 

Combining  equations  (B-l)  and  (A-41)  gives 
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To  a  good  approximation,  equation  (B-2)  can  be  expanded  using  the 
binomial  series  to  give 
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where 


j-<1 

m? 


1  +  u  C. 


( B-3) 


u  i  r 


rj_iV5. 

Ll  *  r  Id  „2. 


(y-1)  m2 
2  Mo 


Specializing  to  n  *  6,  equation  (A-46)  can  be  integrated  using  equation 
(B-3)  to  yield 

V  W’  ]m[l  *  *  "efe)1/6-  ‘P  -]  1  cp -r  ,B-4) 


where 


B3  S 


18g2(6)  (r$  -10'6) 


1  +  r  M2 


1/6 


;  g(6)  -  0.52 


( B-5) 
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The  speed  ratio  ue/u0  can  be  obtained  from  equation  (B-4)  since 

vuo  ■  ^  -  v1/2 

If  sr  denotes  the  point  where  C  *  4/7,  then  from  equation  (B-4) 

C  p 


To  obtain  an  expression  for  C^( s)  for  values  greater  than  4/7,  use  is 
made  of  the  momentum  integral  in  the  streamwise  direction  which  can  be 
written  as 

(0'8) 


where 

9  is  the  momentum  thickness 
H  is  the  boundary  layer  shape  factor* 
rw  is  the  wall  shear  stress 

and  'prime'  denotes  differentiation  with  respect  to  s. 

Setting  the  wall  shear  stress  to  zero  and  assuming  the  shape  factor  is 
constant,  equation  (B-8)  Integrates  to 

/  *  l+H/2 

pe  9  *  °l  /  ^  "  V  (B"9) 

where  a1  is  the  constant  of  integration. 


*  For  turbulent  flow,  separation  is  assumed  to  occur  for  values  of  H 
between  1.3  and  2.4. 


When  Cp  *  4/7  the  inner  region  of  the  boundary  layer  (as  described  in 
Appendix  A  in  the  section  on  the  Inner  Field  Development)  extends  all 
the  way  to  the  edge  of  the  boundary  layer.  In  the  decelerated  flow 
region  where  (ug/uQ  <\j2/7),  we  assume  that  the  speed  ratio  can 
reasonably  be  represented  by  the  first  term  in  the  expansion  given  by 
equation  (A-ll)  evaluated  at  y  »  ye  a  <5.  That  is, 

2 


»)■ 


*  1  •  C, 


4  A2S 


But,  from  equation  (A-38) 


.2  1  pe  1  (dCp\ 


so  that 

■a s~ 


K 

T 


T 


( B-10) 


If  it  Is  assumed  6  remains  proportional  to  the  momentum  thickness  a  for 
s  >_  sc  (a  similarity  solution)  then  equations  (B-9)  and  (B-10)  can  be 
combined  to  give 


ds 


U-cp 

IT" 


„  ,  .  H 

7 


( B-ll ) 


where  02  Is  the  constant  of  propportl  onal  ity  relating  0  to  6  (i.e., 
9  *  o^s).  Integrating  equation  (A-13)  and  selecting  H  »  2  (the  value 
assumed  by  Stratford)  yields 

cp(s)  »  1  -  (a3s  +  a4J"1/2  ;  S  >_  sc  (B-12) 

★ 

where  the  constants  are  determined  by  requiring  that  Cp  and 

★  ★  * 
dC  /ds  (from  equation  (B-12)  match  at  s  *  s-  with  the  C„  and  dC/ds 

p  c  p  p 

expressions  obtained  from  the  separation  criterion.  This  process  leads 
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Momentum  Thickness  Aft  of  a  Normal  Shock 


Consider  the  .development  of  a  boundary  layer  along  a  flat  plate.  Sup¬ 
pose  at  some  point  s$  along  the  plate  a  normal  shock  wave  occurs.  The 
shock  wave  interacts  with  the  boundary  layer  and  modifies  its  velocity 
profile.  As  a  consequence,  the  momentun  thickness  is  substantially 
changed  through  the  shock. 

Figure  C-l  presents  an  illustration  of  a  flat  plate  boundary  layer  in 
the  vicinity  of  a  normal  shock  wave.  Our  immediate  intent  is  to  derive 
an  expression  for  the  momentum  thickness  9  at  station  s  3  ss  +  e  where 
s$  denotes  the  shock  location  and  e  is  a  small  value  measuring  the 
effective  distance  over  which  the  shock  and  turbulent  boundary  layer 
Interact.  The  expression  can  be  developed  by  applying  the  principles  of 
conservation  of  mass  and  momentum  and  the  normal  shock  relations  to  the 
control  volume  shown  in  Figure  C-l.  The  control  volume  is  rectangular 
having  a  length  s  (measured  from  the  plate's  leading  edge)  and  a  height 
<5(s)  equal  to  the  boundary  layer  thickness  at  station  s.  The  conserva¬ 
tion  of  mass  (in  Integral  form)  states: 

<5(s) 

[ 

0 

where  0  is  the  density 

u  is  the  velocity  component  in  the  freestream  direction,  s 
v  Is  the  transverse  velocity  component. 

Quantities  with  subscript  "0"  denote  freestream  values. 


ou]  dy 


3  f  pv d5  +  f 

Vo  Ve 


pvd5 


(C-l) 


Applying  the  integral  form  of  the  conservation  of  momentum  ( s-component) 
to  the  control  volume  gives 
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r«(s) 

J  <p0-Pi  »<Jy  - 


rw(s)<is 
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1 


<S(s) 

{pa2-pQu^)dy  + 


s 

puQvds  + 


( C-2) 


where  Is  pressure  aft  of  the  shock 

tw  Is  the  shear  stress  acting  along  the  wall 
Uj  is  the  s-component  of  the  fluid  speed  aft  of  the  shock  and 
external  to  the  boundary  layer. 


The  pressure  p1  and  the  velocity  component  u^  are  related  to  the  corres¬ 
ponding  quantities  p0,uQ  through  the  normal  shock  relations.  That  is. 


U1  I**1)  Mo  +  2 
uo  ( Y+l) 


(C-3) 


h 

p0 


i  + 


2y 

FFTT 


where  y  •  1.4  is  the  adiabatic  Index  and  MQ  is  the  freestream  Mach 
number. 


Combining  equations  (C-l)  and  (C-2)  gives 

•<S(s)  rs  /-<5(s) 

pvds 


pu(uQ-u)dy 


rs  c 

pvds  *  I 

Je 


(p1-p0)dy  + 


r- 


(s)ds  (C-4) 


To  a  leading  order  approximation,  we  can  estimate  the  momentum  thickness 
aft  of  the  shock  by  assinning  the  Interaction  region  shrinks  to  zero  and 
causes  a  step  change  in  the  momentum  thickness  (see  Equation  (C-7)  for 
definition  of  9)  at  the  shock  location  ss  (see  Figure  C-2).  That  Is, 
substituting  $*ss+e  Into  equation  (C-4)  and  taking  the  limit  as  e+ 0+ 
gives: 

/*4(S+)  r6(  sf)  rS+ 

U°J  PU  1  -  jj-  dy  *  J  (p1-p0)dy  +  J  ys)ds 


(C-5) 
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where  s”  =  Lim  (s  +  e) 
£■*0”  S 


By  combining  equations  (C-5)  through  (C-8),  we  can  obtain  an  expression 
for  the  momentun  thickness  aft  of  the  shock.  If  we  denote  It  by  9$, 
then 


e 

s 


2 

Ty+IT 


(1  -  M 


5+9 
S  0 


(C-9) 


where  6  =  Limit  6(s  +e)  Is  the  boundary-l ayer  thickness  aft  of  the 
5  e+0 

shock.  Equation  (C-9)  states  that  the  jump  In  the  momentum  thickness 
across  the  shock  Is  directly  proportional  to  the  boundary- layer  thick¬ 
ness  aft  of  the  shock  and  to  the  pressure  jump  across  the  shock.  The 
boundary- 1  ayer  thickness  5$  can  be  written  in  terms  of  the  boundary- 
layer  thickness  5Q  before  the  shock  as 


Ss  -.«0  (1+n)  (C“10) 

us 

where  n  is  generally  a  function  of  the  Reynolds  number  Res(=  and 

the  Mach  number  MQ.  Inger  and  Mason  [13]  provide  an  expression  for  the 
jump  in  boundary-layer  thickness  across  a  shock.  For  Reynolds  numbers 
In  the  range  106  ~  10^,  their  development  provides 


n  -  0.8  (Mq-1) 


(C-ll) 


The  boundary-layer  thickness  &Q  can,  In  turn,  be  related  to  the  momentum 
thickness  9Q  by  assuming  a  power  law  expression  for  the  velocity  profile 
and  substituting  It  Into  equation  (C-7).  That  is,  If 


(C-12) 


then  equation  (C-7)  gives 


-84- 


6  ,  (^lHk-t-2)  9  (C-13) 

0  kg 

where  g  Is  a  compressibility  factor  (see  Reference  [8],  Schllchtlng,  pg. 
674)  (g  •  0.713  for  0.8  <  M  <  1.4). 


Substituting  equations  (C-10),  (C-ll)  and  (C-13)  Into  (C-9)  gives 

0  *0(1  +  (1+4M0)(1-M‘2)]  (C-14) 

s  °L  ( y+1)  leg  0  0  J 

This  Is  the  principal  result  of  this  section.  We  have  obtained  an 

expression  for  the  mamentun  thickness  aft  of  a  normal  shock  in  terms  of 

the  known  momentum  thickness  and  flow  conditions  just  before  the 

* 

shock.  For  typical  values  of  M0  ~1.2,  y  and  k,  it  can  be  shown  that 
0-0 

s — 0  ^  Q(1) 


To  apply  this  result  to  an  arbitrary  airfoil  shape,  we  use  the  following 
expressions  (obtained  from  equations  (17)  and  (18)  in  the  text  of  the 


report  for  the  calculation  of  oQ: 


Rfi^H 


1/2 


;  full  laminar  run  case 


( C-15) 


;  turbulent  case 


where  x0  is  the  arc-length  location  along  the  airfoil's  upper  surface 
where  the  peak  velocity  u0  occurs  and  which  is  also  assumed  to  be  the 
location  of  the  shock  wave.  The  reference  condition  [equations  (C-15)] 
therefore  correspond  to  conditions  just  before  the  shock. 
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Equl valent  Flat  Plate  Length  With  a  Normal  Shock 


With  equation  (C-14)  It  Is  now  a  simple  matter  to  obtain  the  equivalent 
flat  plate  length  with  a  shock.  Consider  the  free  stream  conditions  are 
ul»  pl»  ul»  corresponding  to  conditions  just  aft  of  the  normal  shock 
wave.  For  steady  flow  past  a  flat  plate,  the  wall  shear  stress  and  the 
momentum  thickness  are  related  by 


t 

w 


2  de 
-  plul  Is 


(C-16) 


When  the  free  stream  speed  Is  near  the  sonic  speed  (Mi  ~  l)  and  the 
boundary  layer  Is  turbulent,  we  can  use  the  following  empirical  rela¬ 
tionship  [obtained  by  combining  equations  (C-16)  and  (19)]  which  relates 

u.s 

the  shear  stress  to  the  local  Reynolds  number,  Re  *  -±-  : 

s  1 

tw  «  0.036  Re;1/6  ^  Pluj)  (C-17) 

Equating  equations  (C-16)  and  (C-17)  and  solving  for  s  as  a  function  of 
9  gives  the  equivalent  shock  free  flat  plate  length  sso  as 
9  6/5  u  1/5 

sso*  VT0TUZ2t)  (v")  (c'18) 

Since  9$  Is  proportional  to  9Q  and  9Q  depends  Implicitly  on  the  arc- 
length  location  xQ  of  the  shock  along  the  airfoil's  upper  surface  (as 
given  by  equation  (C-15)},  then  equation  (C-18)  provides  the  relation¬ 
ship  between  the  equivalent  flat  plate  length  sso  and  the  shock  location 
xQ  on  the  airfoil.  Figure  C-3  Illustrates  graphically  sso  and  xQ. 
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Inclplent  Separation  With  a  Shock 


The  Stratford-type  separation  criterion  developed  previously  can  now  be 
used  to  calculate  a  zero  shear  stress  pressure  recovery  which  Includes 
the  effect  of  a  normal  shock  on  the  turbulent  boundary  layer.  In  order 
to  obtain  the  pressure  recovery  which  Is  on  the  verge  of  separation,  a 
relationship  Is  needed  for  the  velocity  profile  before  and  after  the 
shock. 

Suppose  the  turbulent  velocity  profiles  can  be  represented  before  and 
after  the  shock  location  as  power  law  expressions.  That  Is, 


k'ihi 


(before  shock) 


(after  shock) 


(C-19) 


(C-20) 


where  n  and  k  are  constants  related  by  the  shock  strength  from  Gadd's 
analysis  [7]  as 


n  -  (k+3) 


- 


(C-21) 


Substituting  the  velocity  profile  (C-20)  Into  the  definition  of  the 
momentun  thickness  equation  (C-7),  gives  [compare  with  equation  (C-13)]: 


6  .  9  .  g  .  0.713 

ng  5 


IC-22) 


Consequently, 


( C-23) 


where  n  =  1/m 


c2  5  (w»T  fZm»l 


The  pressure  recovery  having  zero  surface  stress  can  then  be  obtained 
from  the  separation  criterion  [equation  (A-46)]  by  Integrating  the  equa¬ 
tion 


*  i^n-2) 


*  1/2 


(•*) 


-  g(n)  (Rs110'g) 


sisso’ 


'p  -  n+r 


(C-24) 


•c  *  0.41  Von  Karman's  constant 

k1  *  0.683  Is  an  empirical  constant  obtained  from  turbulent 
flat  plate  experiments 

r  *  0.9  Is  the  recovery  factor  for  a  turbulent  boundary  layer 
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Stratford  Flows  Starting  with  a  Shock 


To  obtain  the  Stratford  flows  when  a  shock  occurs  at  the  location  of  the 
peak  velocity,  the  same  procedure  used  for  the  shockless  case  can  be 

^  —  o 

followed;  namely  to  Integrate  equation  (C-24)  when  0  _<  Cp  <_  and  use 

the  momentun  Integral  approach  forCp>_-j^.  In  procedure  it 

becomes  necessary  to  relate  conditions  at  points  downstream  of  the  shock 
with  conditions  upstream  of  the  shock.  When  this  is  the  case,  one 
assumes  that  along  streamlines  the  Isentroplc  assumption  Is  valid  every¬ 
where  except  across  the  shock.  The  normal  shock  relations  are  used  to 
relate  conditions  along  the  streamline  across  the  shock  location.  Using 
the  fact  that  the  total  enthalpy  remains  Invariant  all  the  way  along  a 

•ft 

streamline  (even  across  the  shock),  leads  to  the  following  Cp  distribu¬ 
tion  [compare  with  equations  (8-4)  and  (8-5)]: 


S  ‘  B  [(4)1/6'l]2/n''  °  -  S  -  T& 


where 


8  2 


r3ng2(n)(R  -lO"6) 


[ft.  *io' 

so 


1/6  2/n 


1+  v*  ^  u2 

1+  r  T  M1 


,  -  Vso 

sso~  V1 


(C-25) 


and 

cj  -  I  -  C«3**  «»rl/z  ;  cjiwr  C C— 26) 

★  A  * 

where  are  determined  by  requiring  Cp  and  Cp  from  equations  (C- 
25)  -  (C-26)  match  at  s  *  ssc.  This  matching  gives  [compare  with  equa¬ 
tion  (5)3: 


2-n 


|jnl\3  8 

r/j*\  .  i 

(!S£\ 

rr)  n 

IUJ  J 

°4  *  (^r)  *  °3  ssc  (C_27) 
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